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O ; RESUMEN 

Encontramos una nueva solucion analftica para las ecuaciones de evolution 
quimica tomando en cuenta la contribution retrasada de todas las estrellas de 
to < 8 Mq(LIMS) como una estrella represcntativa que enriquece al medio 
interestelar. Esta solucion es construida para tasa de formation estelar propor- 

Oi , cional a la masa de gas en un modelo de caja cerrada. Obtenemos incremcntos 

en C/O y N/O cuando O/H aumcnta, comportamiento imposiblc de igualar 

i-Ch ■ con IRA. Nuestros resultados, obtenidos por dos ecuaciones analiticas, son 

muy similares a aqucllos encontrados por modelos numcricos que considcran el 
tiempo de vida de cada estrella. Este modelo retrasado reproduce la evolution 

■+— » i de C/O-O/H y F - O en la vecindad solar. Esta aproximacion analitica es 

C$ \ una herramicnta util para cstudiar la evolution de elementos producidos por 

las LIMS cuando no sc dispone de un codigo de evolution quimica. 

(N 

j> . ABSTRACT 

We find a new analytical solution for the chemical evolution equations, taking 
into account the delayed contribution of all low and intermediate mass stars 



f - ^ . (LIMS) as one representative star that enriches the interstellar medium. This 

solution is built only for star formation rate proportional to the gas mass 

f^i ■ in a closed box model. We obtain increasing C/O and N/O ratios with in- 

OO , creasing O/H, behavior impossible to match with the Instantaneous Recycling 

Approximation (IRA). Our results, obtained by two analytical equations, are 
very similar to those found by numerical models that consider the lifetimes 

K^ , of each star. This delayed model reproduces successfully the evolution of 

5— i ' C/O— O/H and Y — O relations in the solar vicinity. This analytical approx- 

imation is a useful tool to study the chemical evolution of elements produced 
by LIMS when a galactic chemical evolutionary code is not available. 

Key Words: galaxies: abundances — galaxies: evolution — ISM: abundances 

1. INTRODUCTION 

Chemical evolution models are used to describe the temporal variation 
of the gas mass and the abundances of the different chemical elements that 
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are present in the gas. Their importance relies on the fact that the chem- 
ical history of the studied object (i.e. interstellar medium in galaxies and 
intergalactic medium) can be inferred. Moreover it is possible to get chemi- 
cal information about stellar population properties and characteristics of the 
galaxies that we observe in the low redshift Universe. 

The chemical evolution equations, shown by Tinsley (1974) and corrected 
later by Maeder (1992), are relatively complex and can be solved through 
numerical models. There are some analytical approximations that have the 
advantage of predicting the general behavior of chemical elements in a quick 
and easy way but some precision may be lost. 

One of the most well known analytical approximations is the Instanta- 
neous Recycling Approximation (IRA) (Talbot & Arnett 1971) where the star 
lifetimes are negligible compared with the age of the galaxies. This approxi- 
mation has been widely used because it simplifies the solution to the chemical 
evolution equations, however, massive star (MS) lifetimes are on the order of 
10 6 — 10 7 years while the lifetimes for the low and intermediate mass stars 
(LIMS) are on the order of 10 8 — 10 10 years comparable to the lifetime of a 
galaxy. Hence, IRA provides only a very rough approximation for elements 
produced by LIMS. 

Serrano & Peimbcrt (1983) proposed an analytic approximation related 
to the delays in chemical enrichment in N/O-O/H relation assuming N and 
O yields increase with Z . They present closed and open models (with and 
without gas flows, respectively) that takes into account the delay on nitrogen 
production due to LIMS concluding that nitrogen must be an element mainly 
secondary. 

Later on, Pagcl (1989) presented another approximation for open models 
introducing an arbitrary time delay in order to study the chemical evolution of 
element produced by SNIa and LIMS, such as Fe and Ba (through s-process). 
This time delay term makes the stars release the processed material to the 
interstellar medium (ISM) at a single time after the star formation. After 
this single time delay, the contribution of all type of star is instantaneous, like 
IRA. 

The objective of this work is to present an alternative analytical solution 
to the chemical evolution equations that considers the LIMS lifetimes as only 
one group with delay times during the whole evolution. Those delay times are 
different for each chemical clement and are computed based on the character- 
istics of the stellar yields and the stellar population. The delayed contribution 
to the chemical enrichment of elements produced by this type of star at differ- 
ent times should give results with a precision intermediate between the results 
obtained by using IRA and the ones obtained by numerical codes. As an ap- 
plication, this work aims to reproduce the C/O vs O/H and Y(0) histories 
indicated by the stars and HII regions at the solar vicinity. 

This approximation is a simple tool for theoretical and observational as- 
tronomers who need to include chemical aspects in their computations or to 
interpret observational data when they do not have access to a numerical code 
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of chemical evolution of galaxies. Our equation for the mass abundance of el- 
ement i, Xi(t), would replace the popular equation of Z(t) obtained assuming 
IRA in the closed box regimen, erroneously applied for elements produced by 
LIMS (e.g. He, C, N). 

2. THE DELAYED CONTRIBUTION MODEL 

This approximation finds an analytical solution to the chemical evolution 
equations. In the context of a closed box model, we take into account the 
stellar lifetimes (r) of the LIMS as a group where each of its members die at 
the same time. In this approximation, after a star forming burst, all massive 
stars (MS) with initial mass to higher than 8Mq enrich the interstellar medium 
(ISM) instantaneously while every LIMS with m < 8AIq enriches the ISM 
after one single time tl, for the total mass ejected. 

We divide R, the material returned by the stars to the ISM, in two terms: 
Rl and Rm, the returned mass by the LIMS and the MS, respectively: 

Rl = (to — mii)(f>(m)dm 

JMi 

pMsup 

Rm — / (m — mji)4){m)dm 

J8M Q 

where m R is the remnant mass, Mi and M sup are the lowest and highest star 
that contributes to the ISM, and 4> is the initial mass function. 

Taking the galactic chemical evolution equations for the gas mass in a 
closed box model 



dMgas 



-y(t) + E(t) 



dt 

where ^(t) is the star formation rate and E{t) is the gas rate that is returned 
to the ISM from stars that die. In our approximation, we have divided E(t) in 
two parts, the rate where all MS return material to the ISM, regarding that, 
t — t(to) « t: 

pMsup 

E M(t) = (to - m R )^(t - T(m))4>(m)dm = R M ^{t) 

J8M e 

and the rate where the LIMS return material to the ISM, t — r (to) rs t — Tl ,so 
that: 



f8M a 

E L {t)= / (m-m R )ty(t-T(m))(j)(m)dm = R L ^(t-T L ) 

J Ah 

We have considered that the star formation rate is proportional to the gas 
mass with a constant efficiency ^, ^(i) = vM gas (t), therefore, 

^P = -v{l - R M )M gas {t) + uR L M gas {t - tl) 
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The first term of this equation represents the contribution of the MS while 
the second one represents the delayed contribution of the LIMS. 

In order to get the abundance evolution equation, we divide the stellar 
population yields Pi in two parts as we did with R, Pi = Pj Ji + Pui -, with Pl 1 
and P^ defined by: 

f8M 

Phi = I mpi(m)(j)(m)dm 

/>M sup 

Pm( = / mpi(m)(f)(m)dm 

J8M e 

where Pi{m) represents the stellar yield of element i, that is the fraction of 
initial stellar mass processed and ejected by a star of initial mass m. 

Taking the galactic chemical evolution equation for the mass of gas in 
element i, Fi{t) = Xi(t)M gas (t), in a closed box model, 

dF 

^i = -x i (tm)+E i (t) 

Dividing the Ei in two parts as we did before with E(t) we get, 



pM sup 

EMi(t) = / [(m — mji)Xi{t — T(m)) + mpi(m)]'^ (t — T(m))(f)(m)dm 

■>8M Q 

E Mi (t) = [RMX i (t)+P Mi }^(t) 



E Mi (t) = uR M Fi(t) + vP KH M gas (t) 
And for the LIMS: 



f8M 

^Li(t) = / [(m — mji)Xi(t — t(??i)) + mpi(m)]^> (t — T(m))(j)(m)diTi 

J Ah 

E Li (t) = [R L Xi(t - n) + P Li ]V(t - Ti ) 



E Ll (t) = vR L F t {t - n ) + vP L M gas {t - n ) 

We finally obtain: 

dF 

—^ = -v{l - R M )Fi(t) + vP Mt M gas (t) + vR L F % {t - n ) + vP ht M gas (t - n ) 

For each element i, the group of LIMS enrich the ISM at a time Tj after 
it was formed; Ti has to be chosen as the representative ejecting time of 
the population. Equivalently, a representative mass is also required: this is 
defined as the mass up to which the accumulated chemical yield of stars of 
lower mass is equal to the half of the entire yield of the LIMS population. As 
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the stellar yields are different for each chemical element and strongly depend 
on the metallicity, we calculate the value M repr for each chemical element and 
for each metallicity: 

M repr ^ r8Af Q 

mpi(f>(m)dm = — / rrvpi<f>{m)dm 

Ah 2 J Ml 

3. SOLUTION TO THE CHEMICAL EVOLUTION EQUATIONS 

We solve the differential equations for M gas and Fi assuming for the mo- 
ment that the metallicity of the stars that enrich the ISM is constant in the 
evolution, therefore the stellar properties, r, Rm, Pm, Rl,Pl are constant. 

3.1. Gas Mass Evolution Equation 
The equation to be solved is: 

^P = -i/(l - R M )M gas {t) + uR L M gas (t - tl) (1) 

For t < tl, the equation can be solved by IRA, which gives: 

Mga.it) = M sos (0)e-«' < t < t l , 

where £ = u(l — Rm)- 

Having solved for the first interval, < t < tl, we now substitute this 
solution in (1) when tl < t < 2tl obtaining as solution: 

M gas (t) = M ga8 (0)[e-^ + e-«*- Ti )(* - t l )(uR l )] 
Solving for the other intervals: 

iw.)-jw.»-«g '**<'- ye*)'* 

fc=0 



(2) 



for [n — 1)tl < t < utl- 

Taking into account the delayed contribution due to the LIMS, the gas 
mass equation is the product of the term that comes from the contribution 
of the MS (IRA) times by a summation that represents the contribution of 
LIMS that were born at t — Utl and enriched the ISM at the time t. 

3.2. Gas Mass in Element i 
The equation to solve is: 



dF 

—^ = -i/(l - R M )Fi(t) + uP Mi M gas (t) + uR L Fi(t - n) + vP Li M gas {t - n) (3) 

Following a similar procedure to that of previous section: 
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Fi(t) = M ga8 {0)e-^[Xi(0) + vPmA 0<t< 
Solving for the second interval: 



Fi(t) = M ffas (0)[e-« t [X i (0) + vP Mi t] + e-^-^\t - T t )(isR L )[X t (0) + vP M% (t - r ? ) + P L JR L ] 

for Ti < t < 2r,. 

Solving for the other intervals: 



Fi(t) = M gas (0)e-^ ]T ^'^ k ^ )ivRL)]k x(\X t (p) + vP Mi (t - kn)] + kP L jR L ) 

(4) 
for (n — 1)tj < t < riTi- 



3.3. Chemical Abundances by Mass X{(t) 

The evolution equation for Xi(t) is obtained from the ratio of the mass in 
the form of element i and the gas mass: 



M 9as (t) 

For the total gas mass we assume that all the LIMS population ejects the 
material in form of element i at the same time as the material with all the 
elements, i.e., n = tl- Then we finally get Xi(t) dividing eq. (4) by eq. (2) 
where the abundance ratios are independent of the initial value of gas mass, 
M 5OS (0). Therefore 

Xi{t) = El:^(M) (5) 

where G(k,t) = [tizJSZiK^g^ an d Q(k,t) = uP Mi (t - kn) + X,{0) + 
kP Li /R L , valid for t - kr { > 0. 

To solve the approximation of Tj = tl in M gas (t), we take the gas mass as 
the average of the obtained masses of each chemical element: 

M gas (t) = giJ^iM) 

as the gas mass is independent of the element i, and 5 elements are considered 
here (He, C, N, O, and Z). 

We have found that the values of M gas (t) and each M floa (t, r») are quite 
similar in values except for high v and long r,-. In such cases, Xj(t) show an 
artificial secondary raise (see Section 4 for details). 
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3.4. Stellar Properties 

We have chosen an age, an initial mass function, and a yield set in order to 
study the applications and limits of our approximation. We consider 13 Gyr as 
the age of the models, the time elapsed since the beginning of the formation. 
The Initial Mass Function adopted is the one proposed by Kroupa, Tout & 
Gilmore (1993) in the mass interval given by 0.1 < tu/Mq < 80, hence, 
M sup = 80Mq and Mi = lM0,the canonical lowest mass star that enriches 
the ISM in the IRA approximation. 

After trying several stellar yields published in the literature, we selected 
those that provide us with the most complete information in terms of mass, 
mctallicitics and physical parameters such as stellar rotation, mass loss due 
to stellar winds and elements produced during the HBB stage. Once they 
were selected, we calculated the yields of the population based on its initial 
mctallicity. 

For LIMS we use the stellar yields and remnants provided by van den Hoek 
& Groencwegen (1997) that take in account the mass loss due to the stellar 
winds. These yields are characterized by the parameter r\ that represents the 
mass loss, we chose a range for r\ for Z pop = 0.001, 0.004 and 0.02. In order to 
find the delay times of the LIMS group, we have adopted the stellar lifetime 
by Schaller (1992). We have assumed that no LIMS formed with Z pop ~ 0. 

For MS we use the stellar yields and remnants provided by Hirschi (2007) 
for Z pop = 10~ 8 - 0, Meynct & Maeder (2002) for Z pop = 10~ 5 and 0.004, and 
Macder (1992) for Z pop ~ 0.02. We interpolate linearly by mass the stellar 
yields and remnant mass of LIMS and MS, we also extrapolate linearly in 
mass in order to reach m up = 80 Mq. 

We do not interpolate the stellar properties by metallicity and we combine 
MS and LIMS properties according to the initial stellar metallicity. For LIMS 
of Z pop = 10~ 5 we assumed the Z pop = 0.001 yields and remnants. 

Below wc show the Tables 1 and 2 with physical properties derived here for 
our approximation, Rm, Rl, Pau, Phi, and t^. As you can notice, all values 
of Ti arc < 1 Gyr, with the exception of to for Z pop = 0.02, which is 5.89 
Gyr. This difference is important because it causes an artificial behaviour in 
oxygen (only for the case with very high u) and therefore in all abundance 
ratios related to oxygen. 



4. RESULTS 

We study the evolution of mass abundances H, He, C, N, O and Z , using 
the delayed contribution model taking into account primordial abundances, 
X(0) — 0.75, Y(0) = 0.25 for all mctallicitics, and we follow the evolution 
from t = until t = 13.0 Gyr with a time step of At = 0.01 Gyr. In this 
section we do not show the evolution for Z pop ~ 0, since we assume that all 
Pop III stars are MS so that, C/0(t) and N/0(t) are constant according to 
IRA. 
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TABLE 1 



RETURNED MASSES, STELLAR POPULATION YIELDS, 

REPRESENTATIVE MASSES (IN M ) AND DELAY TIMES (IN GYR) 

FOR LOW AND INTERMEDIATE MASS STARS. 



Zjp p 


~ 0.001 


0.004 


0.02 


Rl 


0.217 


0.225 


0.238 


P He x KT 2 


1.261 


1.035 


0.801 


M He 

repr 


2.270 


2.414 


2.646 


T~He 


0.719 


0.635 


0.623 


p c x icr 3 


2.680 


1.717 


0.585 


M c 

■ Lrj ~repr 


2.022 


2.223 


2.761 


TC 


0.968 


0.788 


0.554 


P N x 10" 4 


11.846 


9.300 


8.162 


M N 

repr 


4.769 


5.094 


4.769 


TN 


0.110 


0.099 


0.121 


Po x 1CT 4 


2.937 


1.571 


0.790 


M° 

J.r.i re p r 


2.298 


2.362 


1.185 


TO 


0.697 


0.672 


5.890 


P z x 1(T 3 


4.206 


2.832 


1.492 


M z 

■ Lrj ~repr 


2.546 


2.779 


3.126 


TZ 


0.532 


0.432 


0.390 



Next we discuss the result obtained when using our model for different 
values of the gas consumptions /i = M gas /M t ot as well as for different metal- 
licities. 

4.1. Evolution of C/O vs O/H 

In Fig.l, we show the evolution of C/O vs O/H for three different metal- 
licities and three /i values. In each model the value of fi is reached at t g = 13.0 
Gyr. The horizontal lines represent the results assuming IRA, the lower lines 
when only MS are considered and the higher lines when both MS and LIMS are 
considered. Assuming IRA, the abundance ratios by number can be written 
as, C/0(t)= g$jg and 0/H(t)= 07 ^;„, 

Initially, O/H depends on the star formation rate (SFR), the smaller the 
gas consumption /i the higher the O/H abundance value. 

Massive stars are the first that die and so the C/O values obtained with 
the delayed model are identical to the ones obtained using IRA considering 
only massive stars. The effect of delay can be seen when the curve begins 
to increase its slope making the abundance ratio C/O grow. Afterwards, 
the curve flattens coinciding with the IRA case considering MS and LIMS 
together. 



DELAYED CONTRIBUTION OF LIMS 

TABLE 2 

RETURNED MASSES AND STELLAR POPULATION YIELDS FOR 

MASSIVE STARS. 



ZjpQip 


-0 


1.0 x 10~ 5 


0.004 


Z=0.02 


R M x 1(T 2 


7.345 


7.339 


7.370 


7.438 


P He X 10- 2 


1.208 


1.212 


1.092 


0.929 


p c x icr 3 


2.694 


0.854 


0.820 


2.802 


P N x KT 4 


2.363 


0.429 


1.333 


4.689 


Po x KT 3 


6.261 


6.104 


7.285 


2.953 


Pz x 10~ 2 


0.988 


1.012 


1.154 


0.893 



4.1.1. Case pl = 0.1 

In this case 90 % of the initial gas mass has become stars at the end of 
the evolution and only 10 % remains as ISM (See Fig. 1, Panel a, b, c). 

When the metallicity has values of Z pop = 10~ 5 and 0.004, C/O increases 
for the first time when 12 + log(0/H) ~ 8.0 dex respectively (Panels a and 
b) that means in the times equal to 0.97 Gyr and 0.79 Gyr that are the delay 
times for carbon ejection from LIMS for this metallicities (Sec Table 1). 

When Zp p = 0.02, C/O is higher than the lower Z pop cases because Pc / Po 
are higher for MS at solar metallicity. Moreover when Z pop = 0.02 the values 
for O/H arc smaller than those for other metallicities at the same times. Both 
effects are due to MS of high Z , the metal rich massive star ejects more carbon 
in the stellar wind stage, leaving less carbon behind to be processed in order 
to become oxygen. 

In the other hand, C yields of LIMS are lower at high Z, again for the 
mass loss rate; stellar winds are intense in LIMS and stars have less mass to 
produce heavy elements. This makes less difference between the C/O ratios 
with IRA for MS and MS + LIMS. 

Moreover, for Z pop = 0.02 the first increment of C/O comes when 12+log(0/H) 
~ 7.5 dex (Panel c) corresponding to a time delay for carbon of 0.55 Gyr. The 
curve keeps growing due to the increase of carbon relative to oxygen. After- 
wards, the second increase is artificial (at 12+log(0/H)~ 8.6 dex) and it is 
produced by the oxygen dilution due to LIMS. In the fi = 0.1 case the SFR 
is a quite decreasing exponential function, therefore for t > to the number of 
LIMS that are diluting is higher than the number of MS that are enriching of 
O the ISM . 

4.1.2. Case [x = 0.5 

In this case 50 % of the initial gas mass has become stars at the end of 
the evolution and 50 % remains as ISM, therefore the SFR is lower and flatter 
than in the case /u, = 0.1. As a consequence the values for O/H are smaller in 
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this case than what was previously found at the same times (see Fig 1 , Panels 
d,e,f). In the other hand, the C/O ratios obtained in this case are almost the 
same as for \x = 0.1 because C/O oc Pc/Po, therefore C/O depends on the 
stellar yields not on \x. 

The oxygen dilution effect is present again for Z pop = 0.02 but to a smaller 
degree than for \i = 0.1 due to the SFR behaviour. In the \i = 0.5 fort t> tq 
the SFR is high so there are many MS enriching the gas making difficult the 
oxygen dilution by LIMS. 



4.1.3. Case \i = 0.9 

In this case only 10 % of the initial gas mass has become stars at the end 
of the evolution and 90% remains as ISM, therefore the SFR is even lower and 
practically flat. 

In Panels (g,h,i) the results for this case are shown. This is the lowest SFR 
that we have studied, which implies that O/H ratio is very low too. In this 
case there is no second significant increase in C/O vs O/H, because the SFH 
is the flattest one and therefore the amounts of massive stars and LIMS are 
almost constant during the whole evolution. 



4.2. Evolution of N/O vs O/H 

In Fig. 2 we show the N/O-O/H evolution for the same [i and Z pop values 
as in the previous section. 

The general behavior of N/O ratio is similar to the C/O ratio because 
N and C are produced by MS and LIMS and O is produced by MS mainly. 
Therefore when O/H ratio increases N/O ratio tends to increase when the 
LIMS eject material at the delay time. Since MS produce much less N than 
C, N/O is lower than C/O for early times. Also N/O evolution presents an 
artificial raise due to O dilution, as in the C/O history. 



4.3. Evolution of Helium with Oxygen 

In Fig. 3 we show the helium by mass (Y) vs oxygen also by mass (O) for 
the same /i and Z pop values as in sections 4.1 and 4.2. Here we have assumed 
for all Z pop that the initial Y and O values are 0.25 and 0.00, respectively. 

Since the delay time for LIMS is lower than 0.7 Gyr, AY values are very 
similar to those obtained assuming IRA for MS and LIMS. The slopes of Y(0) 
changes with Z pop due mainly to the Z dependence of oxygen yields through 
the stellar winds. Also Y(0) presents the second significant increase caused 
by the artificial oxygen dilution due the LIMS. 
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Fig. 1. Evolution of C/O abundance ratio versus O/H abundance ratio for metal- 
licities of the stellar population Z pop — 10~ 5 (a, d, g), Z pop = 0.004 (b, e, h), 
Z pop = 0.02 (c, f, i) and three different values of gas consumption fi = 0.1 (a, b, c), 
/i = 0.5 (d, e, f) , IX — 0.9 (g, h, i). The dotted lines represent IRA results when 
only MS (lower lines) or both MS and LIMS (upper lines) are considered. 



4.4. Limitation of the Approximation 

Based on C/O-O/H and N/O-O/H and Y - O evolutions the artificial 
second raise is present only for Z pop = 0.02 and fi = 0.1. This second raise 
is due to a huge delay time and a quite decreasing SFR. According to our 
approximation when we study O(t), LIMS eject oxygen and the rest of the 
elements at the same time, tl = to- Since Pl q is very low, the rest of 
elements produced by LIMS dilute oxygen and there arc not enough MS to 
counteract that dilution. 

This behavior indicates that tl = Tj approximation is not valid for an 
element i mainly produce by MS, when the LIMS delayed contribution of the 
rest of elements is huge compared to the MS contribution of the element i at 
a fixed time. 
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Fig. 2. Evolution of N/O abundance ratio versus O/H abundance ratio for metal- 
licities of the stellar population Z pop = 10~ 5 , 0.004, and 0.02; and three different 
values of gas consumption /j, — 0.1, 0.5, and 0.9, as in Fig 1. 



5. THE SOLAR VICINITY 

The solar vicinity is the place in the Universe with the largest number of 
observations, therefore we use it to test our approximation. However, this 
approximation was developed for a constant metallicity of the stars, therefore 
Pi, Ti and Ri are constant during the gas evolution of a galaxy. This assump- 
tion is not valid for a real galaxy, because the stars in a galaxy form with gas 
metallicity that changes in time. Therefore, the Figs. 1, 2, and 3 are good 
illustrative examples, but not real ones. 

We compute a model where the stars are formed with mctallicities sim- 
ilar to the gas metallicity when the SFR is ongoing. We have made M gas 
continuous and Fi when the gas reach the mctallicities (Z gas ) of the stellar 
population (Z pop ) assumed in this work. 

Formally, when Z^ op < Z gas < Z£+£ we assume that P t = — > — ^ — and 
Ti = ' ' — . Specifically for < Z gas < 10 8 and Z gas > 0.02 we have 
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Fig. 3. Evolution of Helium versus Oxygen, both by mass, for the metallicities 
of stellar population Z pop = 10~ 5 , 0.004 and 0.02 and three different values of gas 
consumption u = 0.1, 0.5 and 0.9, as in Fig 1. The dotted lines represent IRA results, 
when only MS (flatter lines) or both MS and LIMS (steeper lines) are considered. 



assumed the stellar properties for Z 



pop 



1CT 8 and Z 



pop 



0.02, respectively. 

However, for 10 -5 < Z gas < 0.02 we have considered average stellar properties 
between two consecutive Z pop 's. 

With this combination we obtain more realistic results in the sense that we 
are taking into account the evolution of the stellar populations depending on 
its initial metallicity. Another advantage of this combination is that different 
star formation rate efficiencies can be used as a function of time representing 
in a more realistic way the formation of different components of a galaxy. 



5.1. Observational Restrictions 

Now, we study the applicability and limitations of the delayed contribution 
model with different Z pop 's comparing our theoretical results with existing 
observations of the solar vicinity. We define the solar vicinity as the volume 
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contained in a cylinder centered on the Sun with a radius of ~ 1 kpc and 
height enough to reach objects located in the Galactic halo. 

Since C is better known than N, we will test our delayed approximation 
with the C/O-O/H relation in the solar vicinity and our observational con- 
straints for C and O are: 

• H II regions to test the model results at the present time (Esteban et 
al. 2005, Garcia-Rojas ct al. 2004). 

• Halo and disk main sequence stars at different times as past restrictions 
(Akerman et al. 2004). 

• The Sun as a restriction at 4.5 Gyrs ago (Asplund, Grevesse & Sauval 
2004). 

Since the Sun is located at 8 kpc from the Galactic center, we used the two 
HII regions studied by Esteban et al. (2005) and Garcia-Rojas et al. (2004) 
nearest the Sun: Orion Nebula and NGC 3576 at r = 7.46 kpc and r = 8.40 
kpc, respectively. The C/H and O/H values have been increased by 0.10 and 
0.08 dex, respectively, owing to the fraction of C and O embedded in dust 
grains (Esteban et al. 1998). Akerman et al. (2004) show values for 34 F and 
G dwarf stars from the Galaxy halo combining them with similar data of 19 
stars of the disk. 

5.2. Chemical Evolution Models 

We have built four chemical evolution models to reproduce the O/H value 
in NGC3576. Since our approximation was obtained for a closed box model 
with a SFR proportional to M gas , all models presented in this section follow 
those assumptions and moreover consider the same IMF, mass range and 
stellar yields specified in section 3.4. 

Each model is characterized by the approximation used in the lifetime. 
Model 1 is our delay approximation, which assumes that the whole group of 
LIMS is represented by a specific star and its lifetime is considered. Models 
2 and 3 are models in IRA, for MS only and MS and LIMS, respectively. In 
those models no lifetime is considered. Model 4 is obtained using CHEVO 
code (Carigi 1994) that considers the lifetime of each star until leaves the 
main sequence. 

A closed box model with a SFR = vM gas that reproduces a final O value 
has only one free parameter that is the efficiency v. Therefore, the models 1, 
2, 3 and 4 need v =0.13, 0.13, 0.17, and 0.23 in order to get 12 + log(0/H) 
~ 8.82, resulting in values of /j, 0.29, 0.21, 0.21, and 0.51, respectively. 

Since O is produced by MS mainly and Models 1 and 2 assume IRA for 
MS only the SFR is identical for both models, but the value of M gas obtained 
by our delayed approximation is higher than models with IRA due to the 
contribution of the representative mass of the LIMS as a group. The SFR 
for Model 3 is slightly higher due to the dilution of O caused by LIMS at 
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each time and M gaa are nearby identical because the high SFR in model 3 
counteracts the material returned to the ISM by LIMS. In model 4, v and 
/i values are higher than those obtained by the rest of the models because 
CHEVO code needs a delay for each star that forms until its death enriching 
the ISM, therefore the SFR have to be higher than the rest of models to reach 
the same O value. Moreover the M gas values obtained by Model 4 in the last 
9 Gyr are higher due to the delayed contribution of LIMS. For comparisons 
between closed box models with and without IRA see Fig. 12 in Prantzos 
(2007). 

Despite the fact that it is known that the [Fe/H] distribution shown in 
dwarf stars of the solar vicinity is impossible to reproduce with a closed box 
model, we use models 1-4 to check our approximation in C/O-O/H and Y(0) 
relations. 

5.3. Evolution of C/O with O/H 

In Fig. 4, we show the evolution C/O-O/H of the models obtained with 
IRA, CHEVO and our delayed approximation. 

Comparing models with the observational data, we conclude: i) Model 3, 
assuming the canonical IRA for MS and LIMS, produces C/O values that 
are almost constant and much higher than the observed ones for most of the 
evolution, ii) Model 2, assuming IRA only by MS, reproduces very well the 
C/O evolution until 12 + log(0/H) ~ 8 because the ISM in the halo formation 
is determined mainly by MS. Also, this model shows the C/O rise due to the 
high C yield of MS at high Z. However, this increase is not to enough to reach 
the high C/O values present in disk stars and in HII regions, iii) Model 1, 
assuming the delayed contribution of LIMS, reproduces quite well the whole 
C/O-O/H evolution and specifically the increase of C/O abundance ratio at 
high O/H. 

It is notable that these results are very similar to results obtained with 
the numerical code that considers the lifetime of each star (model 4) . 

5.4. Evolution of Helium with Oxygen 

In Fig. 5, we show the evolution of He by mass (Y) and O by mass (O) 
obtained again using the models with IRA, CHEVO, and our delayed approx- 
imation. In order to present the power of our approximation independent of 
initial abundances we have plotted AY and AO. 

Since there is not good observational data for He in dwarf stars in the solar 
vicinity and in the nearest HII regions, we have plotted the Y and O values 
from HII regions in dwarf galaxies and M17, an inner Galactic HII region, as 
representative of the past and future of the solar vicinity, respectively. 

The Y — O relation obtained from models 1 and 4 keep an almost linear 
trend when O < 2.5 x 10 -3 and O < 4.0 x 10 -3 , respectively, with a similar 
AY/ AO. Then the linearity is lost in both models and the slopes are more 
pronounced due mainly to the Z dependence of O yields. The difference in 
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Fig. 4. Evolution of C/O versus O/H using the different models: continuous line: 
the delayed approximation, dashed line: taking into account the lifetime of each 
star, lower dotted line: IRA with MS, upper dotted line: IRA with MS+LIMS, 
Observations: circles: HII regions by Esteban et al. (2005) and Garcfa-Rojas et al. 
(2004), squares: dwarf stars by Akerman et al. (2004), 0: solar values by Asplund 
et al. (2005). 



the changing O is caused by the interpolation assumed between consecutive 
Z p0 p. In CHEVO code a linear interpolation in stellar yields is used when Z gas 
increases with time, while in our delayed approximation we assumed constant 
yields: the average yields between consecutive Z pop . 

Again the delayed approximation is better than IRA and is more similar 
to that obtained by the numerical code. 



6. CONCLUSIONS 



• We have found analytical equations for chemical evolution in the case 
of a closed box model and SFR proportional to the gas mass where the 
delayed enrichment by LIMS is represented by a single type of star. 
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Fig. 5. Evolution of Helium versus Oxygen by mass. Models as Fig. 4. Observational 
data: M17, HII region at r = 6 kpc by Carigi & Peimbert (2008) (filled circle), 
extragalactic HII regions by Peimbert et al. (2007) (filled triangles). 



• The delay of LIMS with respect to the galactic enrichment for He, C and 
N produces an increase on C/O and N/O with O/H and in Y with O in 
good agreement with the results obtained by numerical models. With 
IRA, the C/O and N/O values are constant with O/H in disagreement 
with observed data and with model results that take into account all 
star lifetimes. 



• For \i = 0.1 and Z pop = 0.02, Y(0), C/O, and N/O values show artifi- 
cial secondary raises due to the O dilution produced by LIMS and not 
because of the increase of He, C and N produced by the LIMS. 

• The delayed approximation was probed successfully in the solar vicinity 
reproducing the main trends of C/O-O/H relation shown by dwarf stars 
in agreement with results obtained with a numerical model that consid- 
ers all star lifetimes. That relation cannot be reproduced at all by the 
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instantaneous recycling approximation. 

• The analytical equation (eq. 5) obtained by our approximation is a 
useful tool to know the chemical evolution of those elements produced 
by LIMS when no galactic chemical evolutionary code is available. 

Part of this work was submitted in the Physics Undergraduate Program at 
the Universidad Nacional Autonoma de Mexico. 
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